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EFFECT OF FUSELAGE AND ENGINE NACELLES ON SOME 
AERODYNAMIC PROPERTIES OP AN AIRPLANE WING* 
By Joan Vladea 



With the aid of the method of J. Lotz (reference 1), 
the writer undertook to solve theoretically the lift dis- 
tribution along the span of an airplane wing, when the 
outline of the wing is uneven. This problem arises in 
the case of a mid-wing monoplane with embedded engine na- 
celles. The fuselage and the nacelles were considered as 
aerodynamically profiled, that is, as lift-producing parts. 
The task was therefore to determine not only the disturb- 
ance caused by the fuselage and nacelles, but also their 
share in the total lift of the wing. 

After the above-mentioned calculations, the induced 
velocity, due to the vortex system by which the wing was 
replaced, was also calculated and measured at the points 
corresponding, in an airplane, to the location of the hor- 
izontal tail surfaces. The tail was in fact wanting, so 
that its interference with the fuselage was disregarded. 
In order, however, to allow for the effect of this inter- 
ference, use was made of some of the data published by 
Gorskyi (reference 2). 

As mentioned by J. Lotz, the method discovered by her 
can also be used for wings in which discontinuities occur 
in the profile, angle of attack, or plan contour. The 
funct ions 

f x ( 6 ) =. a g ( g ) s in 6 = Z a gn Bin n 5 
f ? (6) = r° — sin 6 = 2 7 cos 2 6 



*"Uber den Einfluss des Rumpf es |( ( Gondeln) auf einige aero- 
dynamische Eigenschaf t en des Flugels." Z. P. M. , October 28, 
1933, pp. 555-358. 
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fulfill the Dirichlet conditions which are necessary to 
develop a function in the Courier series.. In what fol- 
lows, we undertake to determine the distribution of the 
lift and of the induced drag for the wing represented by- 
figure 1. The profile for |( the wing is the Gottingen 723 
and for the fuselage the Gottingen 350. These profiles 
have c x = 2.64 and ci = 2.53, respectively. The pro- 
file for the region of the nacelles was tested in the 
Aachen wind tunnel. It has ci = 2.64, and the angle of 
attack for c a = 0 is *-\0,5°. If the assembly i.s also 
considered, the following differences between the angles 
of attack are obtained: 
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The coefficients of the Fourier series are determined 
in the usual manner. 

tt _ 6 X 6 s 

a gl = |-j/"fi(6) sin 8dS= ~ / a x sin 2 6d6+/ a a sin 2 8dS + 

83 tt- 5 3 tt 8 a 

+ / a 1 sin 2 8d8+ / a 3 sin 3 SdS+ / aa.sin 3 6d8 + 
62 §3 n-63 

TT-6! TT -1 

+ / a 2 sin 2 8d8+ / cxxsi^SdS; . 

TT-6 2 TT — 5 1 
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By ' iht egrat ion we- obtain: 



a 



i TT 



^ (a x -as) C-(6!-S 2 ) - J-(sin 2 6 x -sin 2 6a) J + 



+ ^ (cci-a 3 ) (6 3 - i sin 2 83). 

TT «5 

We also obtain the coefficients ag^ in the same way. 

a Sn = t? (a * " as) 



|"s in 
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________________ ^ 



The amplitudes 7 3u are determined "by integrating between 
the same values as before. 
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Since the calculation requires a constant Ci over 
the whole wing span, this coefficient was determined as an 
arithmetical mean. 

C1-C13. - ( CO S8l "CO S&2 ) (C11-C3.2) ~ ( Cl 1 - C1.3 ) COS 83. 

The necessary coefficients for the equational system were 
calculated with the aid of these formulas, the unknowns 
being the coefficients cx e n °^ *ke series 



F(S) = V l *d Ci S a en sin n 8 
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Since the convergence of the system was not very pronounced, 
eight harmonics were calculated. : 

For -:the case of .a wing with ..only the fuselage, the so- 
lutions of the. equational system are as follows; 
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With these values in the following points of the span 
6 = 15°, 22^°., 30°, 45°, ,60°, 57-^°, 75°, 80°, 90°, 
the values 

\~ ^en. sin 11 

and 

t 0 n cc en sin n 6 

were determined, which are proportional to the lift and to 
the induced velocity. The sums were multiplied by t 0 /t- L , 
in-:-arde±--to ' ob-tain comparabl e data for fuselages of differ- 
ent lengths. 

It is obvious that the coefficients 0Cen....a-re functions 
of the angles of attack a x and cc 3 . Hence the .distribution 
of the lift and of the induced velocity is also a function 
of the angle of attack. Figure -2 shows that -.the fuselage 
causes a disturbance in the distribution at small angles 
; o-f; -attack,;: that . this- .disturbance becomes., .zero- art,., a certain 
angle; of-, attack., and that, from there on, the . fu.selage 
share is greater than the disturbance produced-, by it. This 
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variat ion th\e 1 if t' distr i.TDutibn . is alsq.,,accompanied by 
a variation in the induced drag; which is proportional to 
the deviation from the' elliptical, distribution. The func- 
tion of the induced "drag does not "become zero .when the 
function of the lift, passes through zero. This interfer- 
ence effect can also "be explained, as follows. It is known 
that ' "•; 

. „ ; p V s " 

• d A = p V r d y = c a -2 — d y . 

that is , 

r = -Z-T *(»)• 

The circulation is therefore proportional to c a *.(£)•• 
Figure 4 represents c a t as a function of cOe for' the 
fuselage and wing. The straight lines intersect at' a 
point where the interference is zero. • The interference 
is proportional to' the difference "between a e for the 
same value of c a t . 

From figure 2 it can then "be seen that the lift' does 
not "become zero in the interval of the fuselage for- zero 
angle of attack. Even in the interval of the fuselage the 
wing induces a lift which results in a positive angle of 
attack. This can happen only when there is an upward cur- 
rent in this region. This is also confirmed "by the dis- 
tribution of the induced velocity. Hence the induced drag 
becomes negative (thrust) in this region and at this an- 
gle of attack (fig. 5). 

In the case of a wing with fuselage end engine nacelles 
the solutions of the equational system are: 
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igure 6 shows that the influence of the nacelles is al- 
bst negligible, due to t±t# :.^hlmpa ; ired'.pro$i;l:e ^ohatp^ct-e: 
is'tics of ■ th e.' w in g even in . tit is r e g i Oai: «'■ ; , T-h e =•_-$ 1 igh %n ess 
Of ■ the 'influence is •■■due -to; the :an'gle of attack -he two-en. 
this portion and the rest of the .wing. ;• ■■ .! 



A -comparison . with- practical :re suits ...was made ;i by - means, 
of wind-tunnel tests. The profile drag was added to the . 
calculated induced drag.-..- This drag was regarded, however, 
as proportional to/the surface area of the model, which 
corresponds to the fact for normal angles of attack, since 
the frictional drag plays an important role at these an-- 
gles (fig. V) . • ' V - : , 

, As a practical result of the. determination of the 
equation for the lift distribution, ' and hence alsb r for the 
free vortices along the 'span of the wing, -the induced ver-> 
tical velocity . in the el evator field was- also calculated',' 
both with and without - the- fus el'eg'e'. A body of rotat ion' ' 
with C-Ottingen profile 36G was- used as the fuselage* 

A few general formulas will now be established. The 
function r = t 0 V c x 2a'en's in 8' will be used ' for the 
circulation distribut ion. ... The point " A (fig. 8) i.s not 
situated on the axis of symmetry of , the, airplane , so. that 
the required formulas have a' general cha'ract er , ' According 
to. ,the 3iot-Savart law., we can write: 
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As customary, we put 
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.. ■ D.i.re,ct. integration of .the expressions, i's- impossible >■' 
"because the ..integrals, are elliptical/. ■ -ti was therefore •'• 
made graphically for the case when ■ ■■ " '■ '•'•"' 

r\ = 1 and'' £ = . 0 . 00 ;■ . 0 . 05; . 0 . 10 ; 0.15;' 0.20; 0.25 

and indeed. for the wing. without fuselage and for the wing 
with, body-of-rotat ion fuselage. The 'accuracy of the draw- 
ing and of the planimetry was. at least 0,001, so that the 
r e suit. s compare very favbrably with the' results of a tedi- 
ous analytical calculation. 

. For comparison with practical results, the angle was 
measured "between the direction of motion at a point where 
the horizontal tail surfaces would pe and a fixed initial 
direction. This angle was measured for the wing without 1 " 
fuselage and with the body of rotation."' The initial direc- 
tion was assumed to coincide with the ' dir ect ion .of the wind 
at zero lift. For the wing with the" G-'ott ingen prof ile 723, 
with which these experiments were likewise performed, this 
zero lift was obtained at an angle of attack of -6°. All 
the other angles were measured from this original direc- 
tion. 



The angles were measured with the aid of a very fine 
silk thread, as shown in figure '9. Silk Was -chosen "be- 
cause of its lightness and "because its su'rfa'ce"' -'is rough 
enough for it to be blown in the direction of the wind. 
The angle, made by the silk thread with the initial direc- 
tion, was measured with a. telescope which was provided 
with a reticle and which could be rotated about its axis. 
The accuracy of this apparatus was 0.02 . In order to de- 
termine the effect of the weight of the thread, the angle 
of the thread was first, measured in the customary suspon- 
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sion of the mo.del in., the wind tunnel and again after turn- 
ing the model 180° about the axis AB. The angular differ- 
differerice thus obtained, -which was • a'bout .'0 . 3- -.'alt the same 
angle of attack, was used in the final determination of 
the angle. The angles were thus measured for the cases 

Ti = 1 and £ = 0.00'; 0.05; 0.10; 0.15; 0.20; 0.25. 

As shown by figures 10 and 11, the experimental re- 
sults agree very well with the theoretical calculations 
far the case .of .a. wing without fuselage. This proves that 
the assumptions were justified. There is no good agree- 
ment in the case of a wing with fuselage, due to the. neg- 
lect- of the lift distribut ion along the length of the- fu- 
selage 'in the theoretical calculations and to the. fail- 
ure to calculate sufficient harmonics. The assumption is 
justified that the whole lift can "be concentrated in the 
centerof pressure,- en which as'sumpti'on the calculation 
of the interference is based. For the case in question,. 
th&; fuselage is much nearer to. the test points, than the 
wing -or the center of pressure". ' Also; the otherwise unim- 
portant- lift along the length of the fuselage has so much 
influence that this assumption is no longer justified. 1 
In the plotted curves, a great reduction in the induced 
velocity in the region covered "by the fuselage can he de- 
tected, the induced velocity being even directed upward at 
small angles of attack. A considerable part of the region 
where these phenomena are noticeable is in the interven- 
ing . space, where the influence of the tail. surfaces is 
very weak (dead-air region). Therefore this disturbance 
phenomenon is not of very great importance. For the re- 
maining region, the induced velocity measured has about 
the same rvalue as for a wing without fuselage. Of course 
these . values are also affected by the shape and width of 
the fuselage. 

^The following formula can be used in the case of a 
wing without fuselage for calculating the induced velo.c-* 
i-ty -along the span of the horizontal tail surfaces, since 
the induced velocity varies b\it slightly in this restrict- 
ed . region. 



1T.A.C.A. Technical Memorandum No. 736 



9 



The first part of the formula yields the induced velocity- 
due to the supporting vortices; the second part, that due 
to the .f_r_.ee vortices. lor the case r\ = 1 1 the part due 
to the supporting vortices is about 1/3 of the tetal in- 
duced velocity and cannot therefore "be disregarded, Of 
course this formula can be used only for fuselages of nor- 
mal shape and width. 

In the above-mentioned work of Gorskyi the shape of 
the curves for the vortices of the induced velocity are 
the same as in our case. He made the measurements with 
the aid of two pitot tubes under 90°, which were therefore 
also very sensitive for the direction of flow. He also 
measured the absolute value of the velocity at the points 
in question, which is not possible by our method. In gen- 
eral he obtained in the region covered by the fuselage a 
stronger upward current than we did. In the remaining 
region, the values are nearly the same. The interference 
between the tail surfaces and fuselage is not at all no- 
ticeable, the results being practically the same for both 
cases with and without tail surfaces. Gorskyi also meas- 
ured the asymmetry in the downward current from the pro- 
peller. 

I am indebted to Dr. C. Wies el sb er ger , the director 
of the aerodynamic institute, for the assignment of the 
subject, for his advice during the work and for the use 
of the wind tunnel for the control tests. 



Translation by Dwight M. Miner, 
National Advisory Committee 
for Aeronautics. 
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Figure 1. -Arrangement of fuselage and engine nacelles. 
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Figure 2. -Lift distribution over span of wing with 
fuselage. 
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Figure 3. -Distribution of induced velocity over span 
of wing with fuselage. 
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Figs. 4,5,6 
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Figure 4.- Circulation for 

wing and fuselage 
plotted against angle of 
attack. 
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Figure 5.- Effective angle 

of attack induced 
on fuselage "by wing at zero 
angle of attack. 
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Figure 6.- Lift distribution over span of wing 
with fuselage and na-celles. 
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Figure 7.-Polars of wing alone and of wing with fuselage and nacelles 
calculated and measured. 
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Figure 9 . -Measurement of angle 
between direction 
Figure 8. -For calculating downward of down current and wing chord. 

current in region of hori- 
zontal tail surfaces. 
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Wing without fuselage 

Figure 10. -Direction of dovm current at horizontal tail surfaces vs. 

span($ of serai span) at various angles of attack of wing. 
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Figure 11. -Direction of down current at horizontal tail surfaces vs. 

span($ of semispan) at various angles of attack of wing 
with fuselage. 



